oo 

o 

O 
(N 



Elasticity and Viscosity of a Lyotropic Chromonic Nematic Studied with Dynamic 

Light Scattering. 
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Using dynamic light scattering, we measure for the first time the temperature-dependent elastic 
moduli and associated orientational viscosity coefficients of the nematic phase in a self-assembled 
lyotropic chromonic liquid crystal. The bend Kz and splay K\ moduli are an order of magnitude 
higher than the twist K2 constant. The ratio KzjK\ shows an anomalous increase with temperature; 
we attribute this to the shortening of the aggregates as temperature increases. The viscosity coeffi- 
cients also show a significant anisotropy, as well as a strong temperature dependence; in particular, 
the bend viscosity is three orders of magnitude smaller than the splay and twist viscosities. 
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Molecular self-assembly in solutions often results in 
anisometric aggregates capable of orientational order. 
The simplest examples are end-to-end "living polymer- 
ization" , formation of worm micelles by surfactants, and 
face-to-face chromonic assembly of disc-like molecules 
into stacks [l|, [2j, La]. In many systems, ranging from 
organic dyes and drugs [l| to DNA [J], this mostly one- 
dimensional aggregation produces a broadly polydisperse 
system of "building units" that form nematic or columnar 
phases and are generally classified as lyotropic chromonic 
liquid crystals (LCLCs). The structural organization and 
properties of LCLCs should depend on the properties 
of both the individual molecules and their aggregates. 
Since the length of the aggregates is not fixed by co- 
valent bonds, it is expected to vary strongly with con- 
centration, temperature, ionic content, etc., making the 
LCLCs very different from the classic liquid crystals in 
which the building units are individual molecules of fixed 
shape. Despite the growing interest in the LCLCs [Hj|, 
very little is known about their molecular structure and 
practically nothing is known about their elastic and vis- 
cous properties. 

In this paper, we use dynamic light scattering from ori- 
entational fluctuations to determine experimentally the 
fundamental elastic and viscous properties of a nematic 
LCLC, following the approach developed by Meyer's 
group for lyotropic polymer nematics [5[. In particu- 
lar, we study 14wt% solutions of Disodium Cromogly- 
cate (DSCG) purchased from Spectrum Chemicals, 98% 
purity, in deionized water (initial resistivity 18 MQcm). 
The homogeneous nematic phase is stable up to the tem- 
perature Tnj « 301 K, above which the mixture is 
in the biphasic nematic - isotropic state. Our experi- 
ment was facilitated by the recent development of align- 
ing techniques for LCLCs [6|] - specifically, the flat glass 
substrates used to sandwich 16^m thick samples were 
coated with a buffed layer of the SE-7511 polymer (Nis- 



san Chemical Inc.) for uniform planar orientation of the 
director n [6( . The cells were carefully sealed with epoxy 
to prevent water evaporation. To mitigate possible ef- 
fects of sample aging, we allowed the cells to equilibrate 
for 48 hours; subsequent measurements were completed 
within 15 hours. Longer equilibration worsens alignment, 
causing parasitic static scattering. The sample tempera- 
ture was controlled within 0.1°C by a hot stage (Instec). 
Light of wavelength A = 532 nm from a diode-pumped 
solid state laser (Coherent, model Verdi V8) was focused 
(beam waist 50/zm) normally onto the sample. The in- 
cident light polarization was vertical; the analyzer direc- 
tion and the scattering plane were both horizontal. 

In the two studied geometries the nematic director n is 
either (1) perpendicular (splay + twist geometry) or (2) 
parallel (bend-twist geometry) to the scattering plane. 
Homodyne cross-correlation functions of the scattered 
light intensity (evenly split between two independent de- 
tectors) were recorded as a function of time on a nanosec- 
ond digital correlator. In geometry (1), two overdamped 
relaxation modes: splay [amplitude A\, relaxation rate 
Li] and twist L4 2 ,r 2 ] are expected [7|. As determined 
by standard selection rules for nematic light scattering, 
the faster mode corresponds to the splay director defor- 
mation, and the slower mode assigns to the twist defor- 
mation. Typical time correlation data are shown on the 
right side in Fig[T] Double exponential fitting of the de- 
cay of bend-twist fluctuations for geometry (2) , displayed 
on the left in Fig[T] combined with an alternative anal- 
ysis by the regularization method (see figure inset) [8(], 
results in an additional two modes with distinctive re- 
laxation rates T3 and T4 and corresponding amplitudes 
A% and At. The best single-exponential fit, shown by 
the thin line in Fig[T] does not match the slower portion 
of the decay. The faster, more intense mode [^3^3] is 
assigned to bend-twist fluctuations, which are expected 
to dominate the depolarized scattering when n lies in the 
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FIG. 1: Correlation functions (open circles) collected at 
T — 294 C for q — 1.64 x 10 6 m _1 in a predominantly bend 
(left) and for q = 1.00 x lO 7 ™,^ 1 in splay + twist geometries 
(right). Double-exponential fits (thick lines) of the correlation 
functions give respectively: amplitude and relaxation rate of 
splay [j4.i,ri], twist L42,r2], bend [A3, Fa], and an additional 
[AtjT^ mode. The thin line for the bend correlation func- 
tion (left) represents a best single-exponential fit, which is 
clearly inadequate. Inset shows the relaxation time spectrum 
obtained by the regularization method 8] for the bend geom- 
etry correlation function; the small secondary peak confirms 
the presence of the additional mode [A4, T4]. 



scattering plane. The weaker mode [A4, T^ cannot be at- 
tributed to a simple leakage of splay or twist fluctuations 
potentially due to slight errors in setting geometry (2) , as 
T.4 substantially exceeds both Ti and T 2 , even when the 
latter two are measured at an order of magnitude larger 
wavenumber q. In this context, we observe from Fig. [5] 
that all detected modes are hydrodynamic (T ~ q 2 ), as 
is normally expected for orientational fluctuations. 

For normally incident light, the amplitude and re- 
laxation rate of the bend-twist mode scale as A 3 ~ 
{K 2 q\ + K 3 qp- 1 and T 3 ~ K 2 q\ + K 3 q 2 [7]. (Here 
|| means along ft.) For small scattering angles, one has 
K 2 q\ <C K 3 q 2 . Assuming for the moment K 2 ~ K 3 , the 

ratio of these two terms is (sin6'b) 2 /4n 2 , which ranges 
from « 0.01 for the lowest scattering angle (6*& = 15°) to 
0.1 for the largest (65°). Here n is the refractive index of 
the solution. In fact, as shown below, K 3 rj (30 — 40)^2, 
and thus the contribution of K 2 q^ to the amplitude and 
relaxation rate of the mode [^.3^3] is utterly negligible 
over the whole angular range studied. Therefore, the dis- 
persion of bend fluctuations can be accurately measured 
while keeping the incident light normal to the sample 
substrates, thus avoiding rotation of the sample in a q- 
scan and potential problems associated with slight non- 
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FIG. 2: Angular dependencies of the rates and ratios of the 
amplitudes of the relaxation modes. Solid lines represent fits 
described in the text. 



uniformity of alignment. 

The variation of the relaxation rates and the amplitude 
ratios with scattering wavenumber q or scattering angle 
9 for the director modes follows the expected behavior 



Ki 



1 K3 tan 
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where r)i with 



(solid lines in Fig. [2]): 

f = ^tan^ h and4 3 - 4Kl tan2 _ ^ 

i = 1,2,3 stand for viscosities corresponding to splay, 
twist and bend deformations, respectively. Since the re- 
fractive index anisotropy An — 0.02 is small [6j, one has 
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-n sin 8b, where 8 st and Q^ are the 



scattering angles, in the splay+twist and bend geome- 
tries, respectively. Good fitting of the relaxation rates 
and the amplitudes indicates that data collected during 
the experiment are not significantly affected by sample 
aging. 

From data for the amplitudes measured on heating at 



fixed scattering angles (8 s t = 40° for splay+twist and 
8b = 15° for bend geometries, respectively) we deduced 
the temperature dependencies of the elastic moduli ratios 
(Fig. [3] a). Absolute values of the moduli (Fig. [3Jb) were 
determined by calibrating the scattered intensity against 
an identical cell filled with the standard thermotropic 
nematic (4'-n-pentyl-4-cyanobiphenyl, 5CB) for a pure 
bend geometry at 294 K. We then used the known Kj, = 
18 pN, An = 0.201 for 5CB H, and the temperature 
dependence of An for our chromonematic [6|, to obtain 
K% = 28 pN at 294 K. Finally we combined the results 
for the elasticities with data for the relaxation rates to 
find the associated orientational viscosities (Fig. [3J3) . 

Within the phenomenological Landau-de Gennes the- 
ory, it is expected that the temperature behavior of 
Ki ~ S* 2 and An ~ S are determined by the temperature 
dependency of the scalar order parameter S pj. In the 
chromonematic case, we find that only the twist mode 

follows this expectation, as the amplitude A 2 ~ £ 
is temperature independent [7j (Fig. [3]c). In contrast, 
i^3/(An) 2 and Ki/(An) 2 both decrease as temperature 
increases. The effect can be qualitatively explained by 
shortening of aggregates with increasing T [10|, which 
should affect K\ and K3, but not significantly K% [111 ]. 
Meyer [3J] predicted that in a self- assembled chromone- 
matic, 
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erage aggregate length [2j 



where L — Lq exp ■ 



is the av- 



2k B T 

Lq is a constant, E a is the 
scission energy needed to split an aggregate into two, ks 
is the Boltzmann constant, and d is the diameter of the 
stack. Qualitatively, Meyer's model describes the data 
well (Fig. [3Jd). It is of interest to explore whether the 
fitting parameters are close to the physical expectations. 

For the analysis, we consider K\/K 2 = -j^r^f-7 (to 
avoid the possible role of S) where the coefficient K\ 
is taken as the value of K 2 /S 2 at T = 294 K with the 
experimentally determined K 2 = 0.65 pN and the typical 
S = 0.7 HiGj. Then with K% = 1.32 pN and d = 1.6 nm 
|13| . we use the experimental data for K\jK 2 to deduce 
L = 0.86 nm, L(T), and E a (as shown in Fig. EldVL 
decreases from 75 nm at 294 K to 34 nm at 300 K ilij . 
which is consistent with the estimate L — 18 nm [10] for 
the isotropic phase at 305K. Furthermore, the scission 
energy shows a weak temperature dependency, changing 
from E a = 3.3 x 10" 20 J at 294K to E a = 2.8 x 10- 20 J 
at 300K, which amounts to (7 — 8)fceT, in agrement with 
estimates for chromonic dyes and DNA oligomers [J, [l2| . 

The case of K3 is the most difficult to discuss as it 
might depend not only on L but also on the persistence 
length P of the aggregates. In the model of rigid rods, 
P — * 00, and K3/K1 cc L/d should decrease with increas- 
ing T if the aggregates become shorter llj , contrary to 
the experimentally observed behavior in our chromone- 
matic (Fig. [3] a). The observed increase of K3/K1 with 
T might be consistent with the model of the flexible rods, 
in which [ll| K3 oc P is independent of L while K\ oc L; 
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FIG. 3: Temperature dependence of the elastic moduli (a, b), 
reduced moduli (c), parameters L and E a (d), and orienta- 
tional viscosities (e), for the studied chromo-nematic. 



however, the data on P are not presently available. 

The viscosity coefficients n sp ; ay , r/twist and r\bend all de- 
crease on heating (Fig. [3^), as in other nematics. It is 
of interest to compare them to the viscosities measured 
(at room temperatures) for the low-molecular weight 
5CB [l5j, and for the lyotropic polymer nematic poly-7- 
benzyl-glutamate (PBG) @ : r] bend (5CB : PBG : DSCG) 
= (0.028 : 0.016 : 0.013) kgm^V 1 and rjtwist (5CB : 
PBG : DSCG) = (0.08 : 3.5 : 20.5) kgrn"^" 1 . Clearly, 
T]bend in all three systems is of the same order, while 
rjtwist in the polymer and chromonematic is two-three 
orders of magnitude higher than in 5CB. The large ratio 
of rjtwist I Vbend ~ 1600 measured for DSCG is not unusual 
from the point of view of the th eory [5| and resemble find- 
ings for the polymer nematics [16(: rjb en d for sufficiently 
extended aggregates is associated mainly with the rela- 
tively easy process of sliding of the aggregates along each 



other, while f]twist should involve rearrangements of the 
aggregates. Interestingly, r)twist / Vbend decreases dramat- 
ically with T, to « 450 at 300K. Such a behavior can 
be tentatively explained by the fact that the aggregates 
become shorter at high T, which should dramatically de- 
crease rjtwist and less so r\bend- It is of interest to explore 
in the future the role of possible breaking and recombina- 
tion of chromonic aggregates in the dynamics of LCLCs. 

An interesting phenomenon observed in the course of 
our experiments is the development at a threshold laser 
power of a strong forward diffraction pattern. Evidently 
a transformation of the sample structure takes place 
above the threshold power, which was about 30 mW at 
T = 294 K, and decreased on heating, vanishing near the 
temperature Tmi = 301 K. To avoid any significant in- 
fluence on our results reported above, we kept the laser 
power as low as feasible between 1.5 and 3 mW and well 
below the threshold for all T studied. This is also why 
we do not include in Fig[3]data measured in the vicinity 
of T/vj, which might be affected by the laser power. 

It remains to discuss the additional mode [Ai,^] ob- 
served in the bend geometry. Since its amplitude is ~ 7 
times smaller than the amplitude of the bend mode, its 
influence, if any, on the results obtained in the bend ge- 
ometry, is small: K3 might be slightly overestimated (by 
less than 15%). Concerning the origin of this mode, we 
can suggest that it corresponds to the translation diffu- 
sion either of monomers, short aggregates (L ~ d) or to 
a transverse sliding of DSCG molecules or their clusters 
within the overall aggregates. From data for I\ shown 
in Fig. [2t> using Einstein's formula for the diffusion coef- 
ficient of spherical scatterers, we find the hydrodynamic 
diameter dh = 3.2 nm, which is twice the aggregate di- 
ameter d « 1.6 nm [13j and could reflect short aggre- 
gates. However the studied system is expected to be 
poly-disperse in length and it is not likely that only the 
diffusion of short aggregates would be detected. On the 
other hand, if we assume that the mode [^4^4] corre- 
sponds to the lateral sliding, we have to account for the 
drag force from the neighboring molecules. According to 
ref. Ljj this can be done by replacing the solvent vis- 
cosity rj with the effective viscosity r] e ff = crj, where the 
drag coefficient c « 2 [l7J. If we use r] e ff instead of r\, 
we arrive at dh ~ d. It is thus plausible that the mode 
[^4^4] corresponds to the sliding of molecules within 
the aggregates. This sliding can provide a bending flexi- 
bility mechanism for the aggregates that differs from the 
known bending mechanisms in other rod- like nematics. 

To conclude, we have measured for the first time the 
temperature dependencies of the elastic moduli and cor- 
responding viscosities for a new type of lyotropic ne- 
matic, the self-assembled chromonic nematic. The ob- 
served anomalous temperature behavior of the splay and 
bend elastic constants indicates that the system is differ- 



ent from the well studied nematic phases of thermotropic 
low-molecular weight liquid crystals and their lyotropic 
counterparts formed either by hard rods or flexible fila- 
ments. One of the key mechanisms contributing to these 
differences is the temperature dependence of the aver- 
age aggregate length. These and other features increas- 
ingly point to the chromoncmatics as a phenomenologi- 
cally richer system than liquid crystals formed by build- 
ing units of a fixed size. 

This work was supported by NSF Materials World 
Network on Lyotropic Chromonic Liquid Crystals, grant 
DMR076290, by Samsung Electronics Corp. and by Fun- 
damental Research State Fund Project UU24/018. 



* Electronic address: odl@lci.kent.edu 
t Electronic address: ssprunt@kent.edu 
[1] J. Lydon, Current Opin. Coll. Int. Sci. 3, 458 (1998). 
[2] P. van der Schoot and M. E. Cates, Langmuir, 10, 670 

(1994). 
[3] R. B. Meyer, in Polymer Liquid Crystals, ed. A. Ciferri, 
W. R. Krigbaum, and R. B. Meyer (New York, Academic 
Press, 1982). 
[4] M. Nakata, G. Zanchetta, B. D. Chapman, C. D. Jones, 
J. O. Cross, R. Pindak, T. Bellini, N. A. Clark, Science, 
318, 1276 (2007). 
[5] V. G. Taratuta, A. J. Hurd, and R. B. Meyer, Phys. Rev. 

Lett. 55, 246 (1985). 
[6] Yu. A. Nastishin, H. Liu, T. Schneider, V. Nazarenko, 
R. Vasyuta, S. V. Shiyanovskii, and O. D. Lavrentovich, 
Phys. Rev. E72, 041711 (2005). 
[7] P. G de Gennes and J. Prost, The Physics of Liquid 

Crystals, 2nd ed. (Oxford, Clarendon Press, 1993). 
[8] S. W. Provencher, Comput. Phys. Commun. 27, 213 & 

229 (1982). 
[9] P. P. Karat, N. V. Madhusudana, Mol. Cryst. Liq. Cryst. 
40, 239 (1977). 
[10] Yu. A. Nastishin, H. Liu, S. V. Shiyanovskii, O. D. 
Lavrentovich, A. F. Kostko, and M. A. Anisimov, Phys. 
Rev. E70, 051706 (2004). 
[11] P. G. de Gennes, in Polymer Liquid Crystals, ed. by A. 
Ciferri, W. R. Krigbaum, and R. B. Meyer (New York, 
Academic Press, 1982). 
[12] V. R. Horowitz, L. A. Janowitz, A. L. Modic, P. A. 
Heiney, and P. J. Collings, Phys. Rev. E72, 041710 
(2005). 
[13] N. H. Hartshorne and G D. Woodard, Mol. Cryst. Liq. 

Cryst., 23, 343 (1973). 
[14] The data remain essentially the same when one uses 
, A \'i for analysis. If one uses K\{T) without S- 
normalization then the fitted L varies from 39 nm to 13 
nm in the same temperature interval. 
[15] M. Cui and J. R. Kelly, Mol. Cryst. Liq. Cryst, 331, 49 

(1999). 
[16] S. D. Lee and R. B. Meyer, Liq. Cryst, 7, 15 (1990). 
[17] G Li and J. X.Tang, Phys. Rev.E69, 061921 (2004). 



